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Abstract 

We study f(R) modifications of Einstein's gravity in a viscous Bianchi 
type I for a Kasner form metric. We get the f(R) function that satisfy 
dominant energy condition. 

1 Introductions 

Observation data[TJ El [3] indicates that there are several shortcoming of 
standard gravity ( SGR)[H [5j [6] related to cosmology, large scale structure 
and quantum field theory. Several attempts preformed to solve these prob- 
lems by adding new and unjustified components in order to give adaptable 
picture with present universe. Almost all of these considerations are mainly 
investigated in a spatially flat homogeneous and isotropic universe described 
by Friedmann-Robertson- Walker (FRW) metric. The theoretical studies and 
experimental data, which support the existence of anisotropic phase, leads 
to consider the model of universe with anisotropic background. Since, the 
universe is almost isotropic at large scale, the studying of the possible effects 
of anisotropic universe in the early time makes the Bianchi I type (BI) model 
as a prime alternative for studding. And solving Einstein equation for a BI 
universe in the presence of dust, stiff matter and cosmological constant done 
in [3 El E] • A kind of approach to be studied in the present work is allow an 
anisotropy in the early universe, and explore where this assumption leads. 
Although, concerning with present universe, it is known that the anisotropy 
is small. Thus Bunnet et al.[10\ have recently done a statistical analysis of 
4-yr data from the COBE satellite, and concluded that the present amount 
of shear (a/H)o is lass than 3 x 10 -9 . In the the early universe, however, the 
anisotropy may have been larger. The possibility of a state of anisotropy at 
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very early times is a very natural to investigate, as an attempt to explain, 
the large local anisotropic that we observe in the universe today in galaxies 
and supergalaxies[llj. Some other recent papers dealing with anisotropic 
cosmologist are [121 [TBI ITU IT5] . 

The present paper organized as follows: In sec. 2, we have reviewed f(R) 
models, the conformal transform between Jordan frame and Einstein frame 
and metric of Kasner form metric. In sec. 3, we have considered the field 
equation of f(R) gravity theory in Kasner form. Sec. 4, is related to con- 
sidering the set of solutions of the field equations. Finally, the latter section 
is devoted to conclusions. 



2 Preliminary 

2.1 /(^Theory 

f(R) theory is a modified theory of gravity, in which the Einstein-Hilbert 
lagrangian i.e. R, is replaced by an arbitrary function of scalar curvature 
i.e. f(R)- So the action for a f(R) theory of gravity takes the following 
form. 



S f(R) = J d^xy/^g + «£m(Vv/>)J . W 

where is the metric of space time, i/j is the matter field and J0. m (g fiu ^) is 
the matter lagrangian. Here we choose k = 1. We assume that the metric 
g^ u is only independent variable and by variation of the action with respect 

V S ives: 
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R^f(R) - ^f{R)g^u - v ll v u f(R)+g^af'(R) = T. 



where TJJ}. 
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2.2 Transformation to a Scalar- Tensor Theory 

As regards, the f(R) theory can be recast as a scalar tensor theory via a 
famous conformal transformation as 



e 



-2/90 



f(R), (3) 



we can rewrite the action (JTJ) as follows 



R 1 

T 



S ST = I d A x^g ( - - -cT^^u4> ~ V{4>) + Cie-W+g^ilH) , (4) 
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where, we define the Einstein frame metric by a conformal transforma- 
tion g^ v = e _2/3<?i <? MI , . Also the potential V(cj)) is given by 

Rf'(R) - f(R) 
V{<i>) = 2ff[R)2 ■ (5) 

Then, variation with respect to the independent variables g^ u and cj> give 
rise to the following field equations 



l -(V<P? + V{<P) 



+t™, (6) 



n<t> = v'(4>) - pT m . (7) 

In the above expression, the covariant derivatives, obey V ^g^v = 0, and 
all indices are lowered and raised with g^ v unless stated other wise. 

2.3 The Kasner Form Metric 

The line element of the Bianchi type-class I can be expressed as follow 

ds 2 = dt 2 - A 2 (t)dx 2 - B 2 (t)dy 2 - C 2 (t)dz 2 , (8) 

where, the metric function, A, B, C, being the function of time, t, only. 
This model is an anisotropic generalization of the Friedmann model with 
Euclidean spatial geometry. Where, the expansion factors, A, B, C, deter- 
mined via Einstein equation. At the moment, we shall restrict ourselves to 
the Kasner form of the metric as: 

ds 2 = dt 2 - t 2pi dx 2 - t 2p2 dy 2 - t 2ps dz 2 , (9) 

where p\, p2, P3 are three parameters which we will require to be constant. 
The space is anisotropic, if at least two of the three Pis(i = 1,2,3) are 
different. 



3 The Field Equation of f(R) Theory of Gravity 
in the Kasner Form Metric 

By making use of ([6]) and (J7|), we can derive the field equation. Considering, 
the Kasner form metric © and introduce the symbols S = p\ + P2 + P3 and 
Q = Pi + P2 the components of Ricci tensor Rqq and Ra obtain as 

follows: 

i?oo = ^ + ^ + 3/30, (10) 
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on the other hands, by considering the energy momentum tensor for a viscous 
fluid as 

T„ v = {p+{P- £0)) u„u v - (P - + 2 V a^, (12) 

where u^, p, P, £ and 7] are the fluid's four velocity, the energy density, the 
isotropic pressure, the bulk and shear viscosities, respectively. The scalar 
expansion and traceless shear tensor are 



(13) 



(14) 



respectively. Here h^ u = u^u u — g^ u is project tensor. The trace of the 
gives : 

R = (Vcj)) 2 + 4V{(/>) - f m , 



substituting, (|T5l) into ([6]) gives 



R 



1. 



So, the components i?oo, Ru follow; 



Ri. 



Roo — e 



2/3<t> 



p-\{ P -Kp-m) 



2 p-(p-tf) + l(f-„ 



23< V((/>) t 2pi 



By using equations (fTU|) . (fTT]) . (fT7|) and (fT8|) one can arrive at. 



(15) 
(16) 

(17) 
(18) 



,2/30 



p--( P -3(p-m 



S-Q Sp<p 



t 2 



+ 



Pi 



+ 



t 2 



+ 



+ 



=2/30 



t 

4t]S 



i + 3/30, 
- 2/3 V + $ 

+ e - 2/3< V(</>). 



(19) 



(20) 
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Since the Pis(i = 1,2,3) are all linearly independent, we get for the shear 
viscosity the following relation: 

^=(I^ + /^) e -W (21) 



Then by substituting, equation ([2T]) into the right hand side of the the ([20]) 

we get 

e 2 ^ ( P - (P - £0)) = 4/3 2 2 - m - - ^1^1 + 2e-**VM, 

(22) 

combining, (|22p and (|19|) . we obtain explicit expression for p and j? as follows 

(23) 

e 2 ^ P = ^(S 2 -Q) + e-%W - 2/3^ + (3/3 2 - i)0 2 . (24) 

Substituting (f23|) and ([24"j) into (J7|), we get the equation of motion for the 
scalar field (j), as follow; 

+ £ _ ^ = _ __ (n0) _ + _ (25) 

where C = (jzf^y (25 - S 2 - Q). 

4 The Solutions 

In this section we want to solve equation (|25p . It is seen that this equation 
has a c/t 2 term (c is a constant), therefore the other terms should phrase 
from order to C(l/t 2 ) term. As a result, we suggest the field <fi to be as 

= 0o Mr-)- ( 26 ) 
to 

By this choosing, we can considering three class solution as follows 
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4.1 Case I 

A class solution of (25) is as 



-2/30 



(1-6/? 2 ) 

7 ,5 



{V {<!>)- W(0)) = o, 

/30 2 - - = 0. 



Solving (|27p . we have 

Then by substituting 
satisfy it, as 



V(0) = y e W . 
in (28), we can obtain an equation which 



{S-l)^> -i3(4 I) ) 2 -C = 0. 
Also we can obtain p, P and rj for this situation as 

Po 



V + 



P 



P 



o 



T T & J- 

t t 2(/34P+l) 



and 



Where 



1-5 



+ 



4P 



t (2/34 7) +l) ' 



(27) 

(28) 
(29) 

(30) 

(31) 
(32) 

(33) 



Pn 



Po 
2 



-S 



S 2 -Q 
2 

Q S 2 



+ 



+ 



(3/3 2 
4 



"4° 



2/35 



1 



3 5/3-2/3-(/3 2 + 2 



4" 



(J) 



2/3^ 



(/) 



3 2 6 

So that by using (j29[) and ([5]) one can obtain the available /(-R) as 

f(R) = -2V + cji?, (34) 

where cj is the integration constant and 2Vq is similar to cosmological con- 
stant A. So that in this class of solution, the function f(R) is similar to 
the standard Einstein general relativity. For the case /3^q ~ 0, which is 
equivalent with f'(R) = 1, this theory is reduce to standard Einstein gen- 
eral relativity theory. In this case the quantity p, P and r] are in agreement 

with \m . 
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4.2 Case II 

The another class of solution is as 



-W4> b 

= to, (35) 



(1 - 6/3 2 ) v ^' H ^" t 2 ' 

t + tj-P^- 1 ^ = 0. (36) 



Where B is the constant, by solving (j35|) we have 



V(t) = V M -k——, (37) 

r 



where 

B(l-6/3 2 )<^ 
2(1 + /3<^ 

Then by substituting (|26p in (|36p . we can obtain an equation which ^q 77 ^ 
satisfy it as 

(S - 1)# 7) - /3(^ /7) ) 2 -(5 + = 0. (39) 
Also we can obtain p, P and r/ for this situation as 

P = — 7rV- + T^- y o, (41) 

r_i + ^(") ) f o M2 (42) 



Here 



Pi = I ^ + ( (3/3 2 - " " 2/35) tf 7) " *) ^° 



Pi = (§5 - | - ^ + (> - 2/3 - (/3 2 + I)tf 7) ) ^ " *) . 
So that by using ([35]) and © we can obtain the available f(R) as 

f(R) = -2V + F R m+1 . (43) 
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Where Fq = w t, 2 \ is the constant and 



m 



V = Vo4 m , (44) 

w = 2S - S 2 - Q. (45) 

For obtain stability of the above f(R), the following condition must be 
satisfied [TT]. 

d 2 f 

d > »• (46) 

for obtaining this constrain and confined the solution power-low in the ex- 
panding universe of Kasner form m should satisfy the following relations [T7] 

m = ^ > 0, w > -2k. (47) 

From (|43|) it is clearly obvious that for ficj)^ 1 = 0, the obtained f(R), (|43|) . 
reduces to f(R) = R which is the standard Einstein general relativity action. 
However, the quantities p, P and rj in these two classes for (3<pQ = are as 

P = ^, (48) 

p _S ±S-S 2 -3Q 

6i 2 ' 1 j 
1-5 

" = -2T- (50) 

The corresponding particular cases are studied in detail in |16j . Considering 
that it dominant energy condition (DEC) imply that — p + j£ < P < p+j£. 
Then, from (|4"8|) . (f49l) and (j50|) . for P < p + | £ we have that 1 - 5 < and 
for of case — p + j£ < P we obtain that —25 — S 2 + 3Q < 0. Thus, from 
these expressions it is easy to see that 77 < and p > 0. 

4.3 Caselll 

In this case we assume B = — C, according to ([39]) (j)Q HI hs as 

,(///) _ 5-1 



' 



(51) 



so fc = _ C(5-l)(l-6^) = _ ft25-g»-Q)(g-l) ^ golutions ^ ??) ^ 

as 

+2(5-2) 

y(t) = y t 4(5 - 1) -fc^ T ), (52) 
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we can obtain physical quantities as 



P = Vo + -^s, (53) 
Pi S 

P = ^ + jZ-Vo, (54) 



2(5-1) 



t 2S- 



(55) 



Here 



= " 2S(S ~ + (3/?2 " i* 1 ^ " fc ) '» (s_1> 

«-(§*-§-* + 5 S < S - D - 2(s - 1) - 0? + - *) <f 

and we can obtain /(i?) as 

/( J R) = -2V + F i? s - 1 . (56) 

Whereas the shear coefficient of viscosity at early time is much bigger than 
the bulk viscosity [IB] , so we take £ = 0. And also according to ([55]) we can 
assume that Vq will be equal to A, therefore we can neglect it with respect 
to other terms in (|53p , (|54p . Hence in this case one can write the barotropic 
equation of state as P = jp. It is well known that the dominant energy 
condition (DEC) imply that 7 lies in the interval [-1, 1]. From (|53p and 
(I54p . and the above assumptions, we can obtain the barotropic equation of 
state as 

_ -fc + §S-3-f + |g(g-l)-2(S-l)-(^ + ^)^ 
7 _ fe+ ^-25(5-l) + (3^-i)^ ' ( j 

The quantity 7 has four independent parameters f3,pis(l = 1, 2, 3). We want 
to obtain the available range of these parameters in which the dominant 
energy condition, — 1 < 7 < 1 and S > 1 is satisfied. We have drown the 
equation of state parameter, 7 versus p% for two set of parameters, a) for 
13 = ^ = 0.41, pi = 0.5 and p 2 = 0.3. b) for /3 = 0.5, pi = 0.5 and p 2 = 0.3 
in figure 1. For these two considerations we obtain two different intervals 
for P3 as (0.424 , 0.69 ) and (0.456 , 0.69 ) respectively. It is seen that 
these two set of parameters and corresponding intervals of p^ are satisfying 
the constraint S > 1 and dominant energy condition. Also, according to[6j 



9 



the second law of thermodynamics imply that ij > hence, considering to 
equation (|55p . it is clearly seen that this condition is satisfied. Therefore we 
derived that this model satisfy the dominant energy condition (DEC) and 
the second law of thermodynamics (SLT), simultaneously. From these two 
sets of data seen that by increasing the value of /?, the constraint S > 1 
, DEC and SLT is satisfied for smaller interval of P3.i.e Apg = 0.266 > 
Apf ] = 0.234. 




0.45 0.50 0.55 0.60 0.65 0.50 0.55 0.60 0.65 

(a) (b) 



Figure 1: Variation the state parameter, 7 with respect to p^: (a); In this 
subfigure we draw 7 for p\ = 0.5, P2 = 0.3 and /3 = ^= = 0.41. We 
have obtained the interval (0.424, 0.69) for ^3 in which the dominant energy 
condition is satisfied. (b);In this subfigure we draw 7 for p\ = 0.5, P2 = 0.3 
and j3 = 0.5. We have obtained the interval (0.456, 0.69) for p% in which the 
dominant energy condition is satisfied. 
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5 Conclusion 



The main purpose of the present work is to study of the f(R) theory of 
gravity in an anisotropic metric, such as Kasner form metric. We have 
assumed that the universe is filled by a viscose cosmic fluid which is endowed 
with a shear viscosity rj and bulk viscosity £. In this work, we have obtained 
the functions of f{R) gravity that satisfied condition stability and expanding 
universe of Kasner form metric. We derived that this model satisfy the 
dominant energy condition (DEC) and the second law of thermodynamics 
(SLT), simultaneously. At last, by plotting the equation of state pamperer, 
7, versus one of the Kasner's index , p%, for two set of parameters, we have 
obtained two intervals for p^, in which the constraint S > 1, DEC and SLT 
are satisfied. Therefore, it is clear that f(R) gravity model as well as satisfied 
solutions, viscous cosmology in scalar theories for Kasner type metrics. 
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